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Abstract 

QED-effects are known to occur in a strong laser pulse interaction with a counterpropagating 
electron beam, among these effects being electron-positron pair creation. We discuss the range 
of laser pulse intensities of J > 5 • 10 22 W/cm 2 combined with electron beam energies of tens 
of GeV. In this regime multiple pairs may be generated from a single beam electron, some of 
the newborn particles being capable of further pair production. Radiation back-reaction prevents 
avalanche development and limits pair creation. The system of integro-differential kinetic equations 
for electrons, positrons and 7-photons is derived and solved numerically. 

PACS numbers: 52.38.-r Laser-plasma interactions, 41.60.-m Radiation by moving charges, 52.38.Ph X-ray, 
gamma-ray, and particle generation 
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I. INTRODUCTION 



The effects of quantum electrodynamics (QED) may occur in a strong laser pulse inter- 
action with a counterpropagating electron beam. In the well-known experiment Q these 
effects were weak and barely observable. If the laser pulse intensity is increased up to 
J > 5 • 10 22 W/cm 2 the QED effects control the laser-beam interaction and result in multi- 
ple pair production from a single beam electron. 

QED strong fields. In QED an electric field, E ^should be treated as strong if it exceeds 
the Schwinger limit: E > Eg = m e c 2 / (\e\Xc) (see [2J). Such field is potentially capable of 
separating a virtual electron-positron pair providing an energy, which exceeds the electron 
rest mass energy, m e c 2 , to a charge, e = — |e|, over an acceleration length as small as the 
Compton wavelength, Ac* = h/(m e c) ~ 3.9 ■ lCr n cm. Typical effects in QED strong fields 
are: electron-positron pair creation from high-energy photons, high-energy photon emission 
from electrons or positrons and the cascade development (see 3J- [4|) resulting from the first 
two processes. 

Less typical is direct pair separation from vacuum. This effect may only be significant if 
the field invariants as defined in [5}], Fi = (B • E), F 2 = E 2 — B 2 , are large enough. Here the 
case of weak field invariants is considered: |F lj2 1 <C Eg, and any corrections of the order of 
Fi^/Eg are neglected (see (J about such corrections). Below, the term 'strong field' is only 
applied to the field experienced by a charged particle. 

QED-strong laser fields. QED-strong fields may be created in the focus of an ultra- 
bright laser. Consider QED-effects in a relativistically strong pulsed field [sj]: 

eA 

|a| > 1, a = -, (1) 

m e c z 

A being the vector potential of the wave. In the laboratory frame of reference the electric 
field is not QED-strong for achieved laser intensities, J ~ 10 22 W/cm 2 jjj], and even for 
the J ~ 10 25 W/cm 2 intensity projected [8|. Moreover, both field invariants vanish for ID 
waves, reducing the probability of direct pair creation from vacuum by virtue of the laser 
field's proximity to ID wave. 

Nonetheless, a counterpropagating particle in a ID wave, a(£), £ = urt — (k-x), may expe- 
rience a QED-strong field, E = \dA/d£\u(£ — pu)/c, because the laser frequency, uj = c/X, 
is Doppler upshifted in the frame of reference co-moving with the electron. Herewith the 
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electron dimensionless energy, S, and its momentum are related to m e c 2 , and m e c correspond- 
ingly, and subscript || herewith denotes the vector projection on the direction of the wave 
propagation. The Lorentz-transformed field exceeds the Schwinger limit, if x ~ E /Es > 1. 
Numerical values of the parameter, x, may be expressed in terms of the local instantaneous 
intensity of the laser wave, J: 

da (S — p\\) 



3 Xc ,c 

x = - 2T (£-pu 



1.4- 10 3 \ 10 23 [W/cm 2 ]' 



(2) 



In the SLAC experiment [l| an electron beam of energy « 46.6 GeV interacted with a 
counterpropagating terawatt laser pulse of intensity J ~ 10 18 W/cm 2 (|a| < 1). A value of 
X ~ 0.4 had been achieved. An increase in the laser field intensity up to ~ 5 • 10 22 W/cm 2 
(|a| ~ 110) with the use of the same electron beam, would allow us to reach a regime of 
multiple pair creation at x ~ 90. 

Radiation back-reaction. The creation of pairs in QED-strong fields is a particular 
form of the radiation losses from charged particles. At high x an intermediate stage in 
the pair creation process is the emanation of a high-energy photon by a charged particle: 
e — > 7, e (in contrast with x < 1 case, in which the "equivalent" photons from the electron 
Coulomb potential mostly contribute to the pair creation - see 4j and the papers cited 
therein). This photon is then absorbed in the strong field, generating an electron-positron 
pair: 7 — > e,p. 

Although the energy-momentum gain from the strong laser field is crucial in the course 
of both emission and pair creation, still a way to quantify the irreversible radiation losses 
may be found. Specifically, in the ID wave field the transfer of energy, AS, from the 
wave to a particle may be interpreted as the absorption of some number of photons, n: 
AS = nhuj/{m e c 2 ). Accordingly, the momentum from the absorbed photons is added to 
the parallel momentum of the particle: Ap» = nhk/(m e c) = nhu/(m e c 2 ). So, both energy 
and parallel momentum are not conserved, however, their difference is: A(S — p\\) = 0. To 
get the Lorentz-invariant formulation, introduce the four-vector of the particle momentum, 
p = (£,p), and the wave four-vector, k = (^,k) for the ID wave field. Their four-dot- 
product, (k ■ p) = u(S — p\\)/c, is conserved in any particle interaction with the ID wave 
field, including its motion, photon emission, pair creation etc. The sum of this quantity, 
J2 {k -Pf)-, over all particles in the final (/) state is equal to that for the particles in the 
initial (i) state: J2 (k ■ Pf) = J2 (k ■ pi). Each term in this conserving sum is positive (we use 



the metric signature (+, — , — , — )). Therefore, any contribution to this sum from a newborn 
particle exacts a contribution from its parent. 

Regarding the high-energy electron beam interaction with the ultra-strong laser pulse, the 
initially high value of x ~ 90 ensures multiple pair creation. The radiation back-reaction, 
however, splits the initially high value of (k ■ p) between all newborn particles. The reduced 
values of (k ■ p) result in smaller values of x oc (fc • p). The cascade terminates, when all 
particles have X < 1 and become incapable of creating new pairs. 

The radiation losses, thereby limit the cascading pair creation. Particularly, emission of 
softer 7 photons even may be described within the radiation force approximation, which is 



traditionally used to account for the radiation back-reaction (see 




The discussed processes are described by the kinetic equations for the involved particles 
(electrons, positrons, 7-photons). For circularly polarized ID wave of constant amplitude, 
the system of three ID integro-differential kinetic equations is reducible to a large system of 
ODEs, which is solved here numerically. 



II. ELECTRON IN QED-STRONG FIELD 

The emission probability in the QED-strong ID wave field may be found in Sections 
40,90,101 in 13], as well as in fl4j . However, to simulate highly dynamical effects in pulsed 
fields, one needs a reformulated emission probability, related to short time intervals (not 
(—00, +00)), which is rederived in Appendix A with careful attention to consistent problem 
formulation. 

Again, the energy, hu', and momentum, fik! , of the emitted photon are normalized to 
m e c 2 and m e c. The four-dot-product, (k ■ p), is the motional invariant for an electron and 
it is also conserved in the process of emission: (k ■ pi) = (k ■ k') + (k ■ pf). A subscript i, f 
denotes the electron in the initial (i) or final (f) state. 

In the ID wave field the emission probability may be conveniently related to the interval 
of the wave phase, which should be taken along the electron trajectory. The interval 
of time, dt, and that of the electron proper time, dr e , are related to as follows: dr e = 
dt/S = dC,/[c(k ■ p)]. The phase volume element for the emitted photon is chosen in the 
form d 2 \i! L d{k ■ k'). The emission probability, dWfi/(d^d(k ■ k')), is integrated over d 2 k.' ± , 
therefore, it is related to the element of the phase volume, d(k ■ k') (see detail in Appendix 
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A): 



dW 



fi 



a (iT K 5/3 (y)dy + KrK 2 / 3 (i 



d(k ■ k')d£ 

(k ■ k') X e 



\/3irX c (k -p., 
(k • k') 



(3) 



Xe 



\{k-p t 



da 



da 



A 



c- 



(k-Pi) ' Xe(k-Pf) 

Here K u (r) is the MacDonald function and a = e 2 /(ch). 

Collision integral. In QED-strong fields we introduce ^-parameter not only for elec- 
trons but also for 7-photons and relate the emission probability to dxy oc d{k ■ k')\ 



X"/ 



l(k-k') 



da 



X 



dW 



c, 



dx-ydi 



Q 



da 



d£ 



rre->-7,e 



j re— >7,e 



^xl 



(4) 
(5) 



K 5 /3(y)dy , 

Here r = x-y/[Xe(Xe — X7)]; Xy < Xe- The electron parameter, X ei is taken for the initial 
state and its value in the final state is Xe — X-y- 

The distribution functions for electrons and photons may be also integrated over and 
k' ± correspondingly. Thus integrated functions are distributed over (k -p), (k ■ k'). We can 
parameterize locally these distributions via Xe oc {k ■ p), Xy oc (k ■ k') and introduce the ID 
distribution functions, f e (Xe) and / 7 (x 7 ). 

The collision integral (see 15|) describes the change in the particle distributions due to 
emission and accounts for the electrons, leaving the given phase volume, dx e , and those 
arriving into it within the interval, <i£ = a\da/ d^\d^ = 2acXedr e / (3Xc)'- 



Sfe(Xe 

d£ 



fe(x)u: 



Xe 

$fy(Xy 
dl 



e— >7,e 

X^X-Xe 



dX ~ fe{Xe) 



Ac 



u xe^x U X, 



v- 



fe(x)U^d X . 



(6) 



Radiation force approximation. One may exclude the emission of softer 7-photons 
with x-y < £ ^ 1 from the collision integral by changing the spans as follows: 



S + fe(Xe 

d£ 



Xe 



UxW^Y-XedX - fe(Xe) J U£%d Xl 



Xe+e 



$ + fy(XyJ 



fe(xW^d X , 



Xy 



Xy > £ - 



The excluded process should be treated separately: 

5 {rf) fe(Xe) d 



da 



dXe 
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Ul f) fe( X e 



(7) 



(8) 



dt J o 

where the expression for the radiation force, 

e 

ul f) = I X,U^d X7} (10) 



o 



is obtained using the standard Fokker-Planck development (see [15|]) of the collision integral 
at small x-y < £: 

e 

J {fe(Xe + X.W^X, ~ feiXeW^) d X , « 


^We( feiXe) lo X " U ^ dX ")- 

The advective term like that in Eq.([H]), once introduced to the kinetic equation, describes the 
electron transport along the characteristic lines, d\ e + d^U^p = 0. This effect is equivalent 
to that from an extra four-force term, (dp* 1 / 'c?r e ) ra d, in the dynamical equation for the electron 
four- momentum, the force being such that: 

Xe ' di dp» \drJ Tad df 
The radiation force is directed along — p^ + j{k ■ p). The two terms describe the energy- 



momentum lost for radiation and those absorbed from t 



re ID wave field in the course of 



emission, their total being orthogonal to p^ (see 11] . 121]). The force magnitude may be 
found from (TTTj) : 



In the first component of this equation the term oc £ controls the radiation energy loss rate, 
Jqed- In dimensional form and related per time interval, /qed = -2am e c 3 U^/(3X c )- At 
Xe < £ 1, /qed tends to the expression for the radiation energy loss rate found within the 
framework of classical electrodynamics. When the radiation force approach is generalized 
for large % e ^> 1, the emission spectrum is modified by the QED effects and only the part of 
this spectrum (which is minor at Xe ^ 1) is embraced by the radiation force approximation. 



III. PHOTON IN QED-STRONG FIELD 



The absorption probability for photons in the ID field is derived in Appendix B. An 
electron-positron pair (e,p) is generated in the photon absorption with the conservation law: 
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{k ■ k') = (k-p e ) + (k-p p ). 

The phase volume element for the created electron, again is chosen in the form d 2 p±d(k-p). 
The absorption probability, dWfi/(d^d(k ■ p e )), is integrated over the transversal momenta 
components and related to the element of the phase volume of electron, d(k ■ p e ), resulting 
in the following collision integral: 

r/e i X£ ' p) = f fMU^'ldx,, (12) 
at; ->xe, P 

5 ~ fM = -/ 7 (X 7 ) r U£%.dx*. (13) 

J 



di 

Here r = x 7 /[Xe(X7 - Xe)], Xe = X 7 ~ Xv < X 7 and 

V3 



£T 7 -s-e,p 



K 5 /3(y)dy 



(14) 



IV. SOLUTION FOR KINETIC EQUATIONS 



As long as the distribution functions are integrated over the transversal components of 
momentum and expressed in terms of the motional integrals, (k ■ p e , p ), their evolution is 
controlled by the collision integrals: 

dfcp^i^ (k • Pe,p,y)) 



{6^ + 5 + + 6-)f eiP!T (15) 



OH 

The derivatives, <9/<9£, are taken at constant (k-p). Eqs. fTT5|) are easy-to-solve for the ID wave 
field of any shape, however, for circularly polarized wave of constant amplitude the solution 
is especially simple. In this case (k-p) are different from x by a constant factor, and Eqs. f|T5|) 
may be solved with derivatives, d/d£, at constant x f° r the functions, f e , P ,-y(£,,Xe, P ,-y)- 

We solve Eqs. (I15j) numerically, by discretizing them at a uniform grid, x = -^^X? I = 
1, 2, 3..., N, with the choice of Ax = 0.1, e = Ax/2. The ^-dependent distribution functions 
at this grid obey the system of 3iV ODEs, which is integrated numerically. At initialization, 
electrons with f e = 5(x e — Xo) ? Xo = 90, counterpropagate in the circularly polarized wave 
field with \d&/d!; \ = 110. This choice corresponds to the SLAC electron beam and the laser 
intensity of J ~ 5 ■ 10 22 W/cm 2 for A = 0.8/im, to be achieved soon. 

In Figfl] the beam-wave interaction is traced during ^- = 5 cycles of the incident laser 
pulse (~ 13 fs). The initial beam electron energy is rapidly converted into 7-photons with 
high x 7 , which then rapidly produce pairs, the typical rates of the processes being of the 




FIG. 1: Distribution functions of electrons and positrons, f e ,p(x)i an d a spectrum of emission, 
Xyfi(x)/Xo, after the interaction of 46.6 GeV electrons with one cycle (left panel) and five cycles 
(right panel) of a laser pulse of intensity J « 5 • 10 22 W/cm 2 (so that % ~ 2£[GeV] — see Eq.fl2J)). 
Here f e — f p is the distribution of the beam electrons and / (f e — f p )dx = 1- 




e/27i 

FIG. 2: Pair production for longer pulse durations, measured in cycles. Other parameters are the 
same as in Fig(TJ 

order of the inverse light period. However, the larger fraction of the new particles is born 
at X < 1, with strongly reduced pair production rate. Slow absorption of photons with 
X 7 ~ 1 — 2 maintains pair production even after tens of wave periods, as shown in Figj2j 

V. CONCLUSION 

We see that the laser-beam interaction may be accompanied by multiple pair production. 
The initial energy of a beam electron is efficiently spent for creating pairs with significantly 
lower energies as well as softer 7-photons. This effect may be used for producing a pair 
plasma. It could also be employed to deactivation after-use electron beams, reducing radia- 
tion hazard. 

The way to solve the kinetic equations is accurate and it does not employ the Monte-Carlo 



8 



method. The solution can be used to benchmark numerical methods designed to simulate 
processes in QED-strong laser fields. 

We acknowledge help and advice we received from S. S. Bulanov, M. Hegelich, 
J. G. Kirk, H. Ruhl, T. Schlegel and T. Tajima. One of us (I.S.) is supported by the 
DOE NNSA under the Predictive Science Academic Alliances Program by grant DE-FC52- 
08NA28616. 



VI. APPENDIX A. ELECTRON IN THE QED-STRONG FIELD: EMISSION 
PROBABILITY 



In weaker fields, especially for the particular case of a harmonic wave, the emitted power is 
given by an integral over many periods of the wave. In this case, the solution of the emission 
problem in the weak harmonic wave field (under the requirement on the wave amplitude 



opposite to that formulated in Ineq.flTJ)) is given in Section 101 of [13| as a sum over multi- 
photon-orders, resulting from the Fourier-series expansion for the (infinitely long) periodic 
wave. This standard approach, however, may become meaningless as applied to ultra-strong 
laser pulses, for many reasons. These pulses may be so short that they cannot be thought 
of as harmonic waves. Their fields may be strong enough to force an electron to expend its 
energy on radiation in less than a single wave period. However, an even more important 
point is that the radiation loss rate and even the spectrum of radiation is no longer an 
integral characteristic of the particle motion through a number of wave periods: a local 
dependence of emission on both particle and field characteristics is typical for the strong 



fields. The latter statement may be found in 



131 ] . Section 101. For the particular case of the 



ID wave field the evaluation of the formation time for emission is provided in [12[ within the 
framework of classical electrodynamics. It is shown that the formation time is much shorter 
than the wave period as long as Ineq.((T]) is fulfilled. 

Here the emission problem is discused for QED-strong fields. We consider a ID wave field 



taken in the Lorenz gauge 



16|: 



^ = ^(0, £={k-x), {k-a) = 0, 

a M = (0, a), and x^ being the 4- vectors of the potential, the wave and the coordinates. 
Herewith the 4-dot-product is introduced in a usual manner: (k ■ x) = k^x^ = ut — (k • x) 
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etc. Space- like 3- vectors (i.e., the first to the third components of a 4- vector) in contrast 
with 4-vectors are denoted in bold, 4-indices are denoted with Greek letters. Recall, that a 
metric signature (+, — , — , — ) is used, therefore, for space-like vectors the 3D scalar product 
and 4-dot-product have opposite signs, particularly: 

/ da\ 2 ( da s 2 



Ml \H - a 



A. Transformed space-time 



A method facilitating many derivations involves the introduction of a specific time-space 
coordinate frame. Introduce a Transformed Space-Time (TST) : 

x ' 1 = ( c t =fx||)/V2, x 2 ' 3 = x ± , 

subscript _L denoting the vector components orthogonal to k. The properties of the TST 
provide a convenient description for the classical motion of an electron in the ID wave field. 
Note first that 

o Adf o _ K k 'V) , j \ £~P\\ 



dx" = p" = v , (p ■ k) 



y/2' r y/2 ' A 

Second, the generalized momentum components, p° and p_i_ = P± + a , are conserved. Third, 

the metric tensor in the TST is: 

G 01 = QlO = i q22 = G 33 = _i Qp, = Q 



[IV 



Note the unusual off-diagonal structure of the metric tensor, resulting in a strange rela- 
tionship between contravariant and covariant coordinates: x° = x±, x 1 = Xq. Specifically, 
the (contravariant) component of the electron momentum, p°, is a motional invariant, as 
long as the vector-potential (and the Hamiltonian, H) does not depend on the (covariant) 

coordinate, x (despite its dependence on x°): 

dp dH 
dt dx 

Finally, the identity, S 2 = p 2 + 1, being expanded in the TST metric, gives: 

1 = 1 + pj = 1 + (p±o - a) 2 
P 2p° V2X(k-p) 
The derivative over x° or, the same, over £ is conveniently related to the derivative over the 
proper time for electron: 



dr 



- + ( — 



c{k-p)± (16) 
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B. Classical trajectory and momenta retarded product 



Many characteristics of emission may be expressed in terms of the relationship between 
the 4-momenta of the electron at different instants: 

^=r'(6)-^ + 2 ^y a) V , (17) 

where 

As a consequence from Eq. ffTTl) . one can obtain the expression for the Momenta Retarded 
Product (MRP): 



W f ) . p( f l)) = i_(^! = 1 + M, (18) 



Note, that the MRP is given by Eq. (fT8]) for an arbitrary difference between £ and £1, 
but only for the particular case of the ID wave field. However the limit of this formula as 
|£ — £i| — > 0, which is as follows: 

2 



da 



or, in terms of the MRP in the proper time, r: 



(P(r)-p(r + 5r)) = l-(5rf { l-^, (19) 

has a much wider range of applicability. Eq.( fl9]) is derived from the equation of motion: 

dp 11 f 

dr m e c 

using the identities: 

(p(r)-p(r)) = l, (p(r) ■ f(r)) = 0, 

d(or) \dr dr J 

Here / M is the Lorentz four-force: 

//* = (/°,f( 3 )) J /° = eE-p, f (3) = eSE + ep x B. 
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C. A solution of the Dirac equation 



The Dirac equation which determines the evolution of the wave function, tjj, for a non- 
emitting electron in the external field, reads: 



iX c ( 7 .|L)-( 7 .a) 



(20) 



7 M being the Dirac 4*4 matrices, (7 , 7 1 , 7 2 , 7 s ). The relativistic dot-product of the Dirac 
matrices by 4- vector, such as (7-0), is the linear combination of the Dirac matrices: (7-0) = 



„o„o 



7 a — 70 — 7 a 



„2„2 „ ; 3„3 



7 . Such a linear combination, which is also a 4 * 4 matrix, may be 



multiplied by another matrix of this kind or by 4-component bi-spinor, such as i/j, following 
matrix multiplication rules. For example, (7 • a)ip is a bi-spinor, as is the matrix, (7 • a) 
multiplied from the right hand side by the bi-spinor, ip. 

The solution of Eq. fpZO"]) in a form of a plane electron wave can be conveniently expressed 
in terms of the classical solution: 



N 



u(p(0) exp(-z / ((P(0 + eA) ■ dx)), 



(21) 



the normalization coefficient N = const. By expanding the phase multiplier in the TST, a 
more convenient form can be provided: 



%p = = exp 



N 



,P_lo • x_U 



X(k-p)x 1 

~7T~ 



X 



c 



(22) 



Here it(p(£)) is plane wave bi-spinor amplitude, which satisfies the system of four linear 
algebraic equations: 

(7'P(0)«(P(0)=«(P(0)> ( 23 ) 



as well as the normalization condition: uu = 2. The ^-dependent phase multiplier, P(£), is 
as follows: 

n 1 + p±(6) 

or, 

<2 ■ (24) 



P®=*p(-±f 



Using Eq. (11711 . one can find: 



2(k-p) di ' 

c i + pi(6 



2 . 



u(p(0) 



1 + 



2(fc • p) 
(7 • k) ( 7 .[ Q (0-a(6)]) 



2(fc.p) 



«(p(eo) 



(25) 
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and verify that Eq. ([22]) satisfies the Dirac equation. To prove all these assertions, note that 
the plane wave bi-spinor amplitude once expressed in terms of the classical solution Eq.flT] 
satisfies the following commutation rule: 



Kt-KO)±i] 



i 



2(k-p) 



(7 ■ 5a) 



2(k-p) 



(7 ■ 5a) 



K7-p(ei))±i] 



(26) 



The latter can be proved using the commutation rules as in Eqs.(22.5) from [13j for the 
Dirac matrices as well as the identities, (k ■ A) = (k ■ k) = 0: 



(t-p(6))-(t-^) + 



2(p(6) • Sa) - (SaY' 
2(k-p) 



(7 ■ k) ± 1 



(7-p(£0)-(7'<fa) + 



2(p(ei)-^)-(7-p(6))(7-^) 



1 + 



2(fc • p) 
(7 - A) (7- 5a) 



(7 -ft;) (7 -fa) 
2(Jfc-p) 

(7-fc)(7-(^a) 
2(fc • p) 



[(7"P(6))±1] 



2(fc-p) 

Thus, Eq.( |26l) is proved. This commutation rule allows us to verify that u(p(£)) is indeed a 
bi-spinor amplitude of a planar wave with the four- vector satisfying Eq. (l23l) . Finally, to 
verify that Eq.(l22l) gives the solution of the Dirac equation, one should account for Eq.(l23l) 
as well as the identity, 



(7.*) 



(7 • fc)(7 • k) 7 



u(p ) = 0, 



which is valid, because (7 ■ k){^ ■ k) = k 2 = 0, 

The advantage of the approach used here as compared to the known Volkov solution (see 



Section 40 in 



131 ]) is that the wave function in Eqs. 



is described in a self-contained 



manner within some finite time interval, (in fact, this interval is assumed to be very 

short below) in terms of the local parameters of the classical trajectory of electrons. This 
approach is better applicable to strong fields, in which the time interval between subsequent 
emission occurrences, which destroys the unperturbed wave function, becomes very short. 



D. The matrix element for emission 



The emission problem is formulated in the following way. The electron motion in the 
strong field may be thought of as a sequence of short intervals. Within each of these intervals 
the electron follows a piece of a classical trajectory, as in Eq.l fPTt) . and its wave function (an 
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ct 



FIG. 3: The volume over which to integrate the matrix element while finding the emission probabil- 
ity: in the standard scheme for the dipole emission (in dashed lines) and in the TST (in solid lines). 
Arrows show the direction, along which the Heisenberg operator advances the wave functions. 

electron state) is given by Eq. (!22|) . The transition from one piece of the classical trajectory 
to another, or, the same, from one electron state to another occurs in a probabilistic manner. 
The probability of this transition, which is accompanied by a photon emission is calculated 
below using the QED perturbation theory. We consider the laser wave field as a classical field 
(Ac(/c • p) < 1) and apply first-order perturbation theory with respect to the Hamiltonian 
of the electron-photon interaction, H oc (J • A). 

The only difficulty specific to strong pulsed fields is that the short piece of the electron 
trajectory is strictly bounded in space and in time, while the QED invariant perturbation 
theory is based on the 'matrix element', which is the integral over infinite 4- volume. 

To avoid this difficulty the following method is suggested, which is analogous to the dipole 
emission theory as applied in TST. Introduce domain, A 4 x = (Ax 1 * S±) * Ax°, bounded 
by two hypersurfaces, £ = £_ and £ = £ + (see Figj3]). The difference £ + — is bounded as 
described below, so that A 4 x covers only a minor part of the pulse. A volume, 




is a section of A 4 x subtended by a line t = const. 

With the following choice for the normalization coefficient in Eq.( l22|) : 




the integral of the electron density in the volume V 
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is set to unity, i.e. there is a single electron in the volume V. This statement follows from 
Eq. (l22p and the known property of normalized bi-spinor amplitudes: u ■ 7 ■ u = 28. Here 
the hat means the Dirac conjugation. 

For a photon of wave vector, (k')^, and polarization vector, introduce the wave func- 
tion: 

, _ exp[-i(k' •x)/X c ] 



(A' 



or, by expanding this in the TST: 



i(k' ± ■ x ± ) iX(k • k^x 1 



A 



where: 



Pp(£) = ex P 



c 



V2X 



c 



1". 



rl ^2 



'2{k-k')X c _ 

Here the photon momentum and photon energy are related to m e c and m e c 2 correspondingly, 
or, equivalently, dimensionless {k'Y equals dimensional {k'Y multiplied by Xc- The choice 
of the normalization coefficient, 



u'V 



p 2nhcXc 

corresponds to a single photon in the volume, V. 

The emission probability, dW, is given by an integral over A A x: 



dW 



aLfLp 



he 



V>/(7 • (A')*)^ i dx°dx 1 dx 2 dx 3 



Here 



Vd 3 k' hcN p d 2 k' ± d(k ■ k') 



(2nX c ) 3 (2nX c ) 2 (k-k>) 



(27) 



(28) 



) 3 

is the number of states for the emitted photon. The transformation of the phase volume as 
in Eq. fTSBI) is based on the following Jacobian: 



OK 



to 



(k> ■ k) ' 



y(9(/c /c) J j c '^ = cons t 

which is also used below in many places. A subscript i, f denotes the electron in the initial 
(i) or final (/) state. The number of electron states in the presence of the wave field, Li f, 
should be integrated over the volume V 



ij 



(2ttX c ) 3 Jv 



PlJ 2(2nrX 3 c (k-p) u 
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E. Conservation laws 



The integration by dx 1 dx 2 dx 3 = c\/2dtd 2 x± results in three S— functions, expressing the 
conservation of totals of and (k ■ p), for particles in initial and final states: 

P±i = P±f + k'_L, (k ■ = (k ■ p f ) + (k ■ k'). 

Twice integrated with respect to dx , the probability dW is proportional to a long time 
interval, At = Ax 1 /(cv / 2), if the boundary condition for the electron wave at £ = is 
maintained within that long time. On transforming the integral over dx° to that over g?£, 
one can find: 



...d 4 x 



(27rX c ) 3 S ± cAtX 



x 



x5 2 ( P ± 4 - Px/ - KMk ■ Pi ) - (k ■ p f ) - (k ■ k')). 

To take the large value of At seems to be the only way to calculate the integral, however, the 
emission probability calculated in this way relates to multiple electrons in the initial state, 
each of them locating between the wave fronts £ = and £ = £ + during much shorter time, 



5t(Z-,Z+) = {l/c) J S t (Od&/(k- Pt ). 
For a single electron the emission probability becomes: 



(29) 



Using 5— functions it is easy to integrate Eq. (l2"T]) over dp±fd(k ■ pf] 



d(k ■ k')d 2 k' ± ~ (47rX c ) 2 (k ■ k')(k ■ Pi)(k ■ p f ) 



where 



T(£) 



cxp 



^c(k-Pf,i) 



Pi(£) and Pf(0 are the electron phase multipliers, P(£), for the electron in initial and final 
states and 



-Jfe". 



(30) 



If in the expression for T(£) the numerator of the fraction is taken for the initial state of 
electron, then the denominator should be taken for the final state and vise versa. 



1(3 



Prior to discussing Eq. (I30j) . return to Eq.( lT7|) and analyze it component-by-component 
in the TST. It appears that three of the four components of that equation describe the 
conservation of (k ■ p) and pj_ = p_i_ + a for electron in the course of its emission-free 
motion. At the same time, yet another component of Eq.( lT7|) . specifically, p l , directed 
along k^, describes the energy-momentum exchange between the electron and the ID wave 
field, maintaining the identity, (p ■ p) = 1. Now turn to Eq. (l30|) . Again, three of the 
four components express the conservation of the same variables in the course of the photon 
emission, while the p l component, directed along k^ describes the absorption of energy and 
momentum from the wave field in the course of the photon emission. Note, that in the case 
of a strong field, the energy absorbed from field is not an integer number of quanta, and 
that for short non-harmonic field it is not even a constant, but a function of the local field. 



F. Calculation of the matrix element: a product of the Dirac matrices 



To calculate the matrix element, one can re- write it as the double integral over d£d£i. 
The integrand in this double integral includes the phase multiplier and the following matrix 
product: 

fi(p/(£))(7 • i>(R(0) [a(p/(6))(7 • i>(ft(£i))r = - (3i) 

Using Eq. (j251) one can reduce the multipliers of the kind of u(pi(£)) <g> u(pi(£i)) to 
u(pi(£i)) <8> u(pi(C,i)), which can be expressed in terms of the polarization density matrix 
at the position £1. Although in a strong wave electrons may be polarized (see 171]). in the 
present work the emission probability is assumed to be averaged over the electron initial and 
final polarizations. Therefore, the polarization matrix is used in the form of 1/2 [(7 •£>,;(£) ) + l] 
and one can re-write the integrand as follows: 



4 



x[( 7 -p/(0) + l] (7-H>< 



Sp{ 



1 + 



■ (7-fc) 
2(fc. P/ ) 

(7-*) 



X 



2(k- Pi 



(7 • HO - a(6)D 



[(7-^(6)) + !] (7-0) 



In transforming the trace of matrix in addition to the mentioned above relationships and the 
Dirac matrix algebra we use the commutation rule Eq.f |26|) . the conservation law Eq. (l30|) . 
and the identity (/*•/) = —1 for a spatial unity vector of the emitted photon polarization, 
which can be also written as (7 • /*)(7 • Pi,/) (7 • = 2(Z* • Pi,/)(7 • I) + (7 • Pij)- Move 
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(7 * P/(0 + 1) to the first position using the commutation rule Eq. (126]) and then move it to 
the last position. In the product [(7 • Pi(£i)) + 1](7 ■ 0[(7 ' P/(£i)) + 1] we should keep only 
odd powers in 7, as long as in the preceding multiplier only odd powers of 7 are present: 



fSp{ 



(7-jO 
2(k-p f ) 



(7 " HCi) ~ a(01) 



( T -r)x 



(7 • k) 



(7" W0-a(6)]) 



[(7-R(6))(7-0(7-P/(6)) + (7-0]} 



2(A;.pi) 

In the product of the first three multipliers separate the even powers in a, which are (7 • 
/*) - (7 • k)(l* ■ k)(a(£) - a(£i)) 2 /[2(/c -pi){k ■ p/)]. The traces of the products of two or four 
Dirac matrices are calculated using Eqs. (22.9-10) from 131 ] . 

(r ■ k)(l ■ k) 



1 + 



(«(0 - 0(6))' 



2{k- Vi ){k-p f ) 

(i* ■ k)(i ■ ft (eo) , (/••*)(/• p/ceo) 



(a(0 - a(£i)) 2 + 



2(fc-pi) 2(fc- P/ ) 
In one of the two terms linear in a we move (7 • I*) from the first to the last position: 

+ 8(F^) Sp {(7 ' k)il ' [a(0 " a(6)])(7 ' ft &M 2 ( r ' p/&))(7 • + (7 

"s(^) Sp {(7 ' fc)(7 ' [a(0 " a( ^ i)])[2(r • p ^ i))(7 • z) + (7 • p*(^))](7 • p/(eo)} 

Calculate traces for the products of four Dirac matrices in the terms linear in a: 

{I* ■ k)(l-k) 



= [(Pi(fc)-P/(fc))-l] 



2(k- Pi )(k-p f ) 



(i*-*)(i-Pi(ei)) + (i*-fc)(i-P/(ei)) 



2(A;-p i 



2(*.p/) 



a(0 - «(£i)) 2 + 



+{-2(K0-«(6)]-p,(6))+ 



(Mei)+Mei)]-n(^-KO-a(ei)]) + 



([a(e)-a(ei)]-ft(ei))(*-o(^-p/(ei)) 



aa(e)-a(ei)]-p/(ei))(fc-o(i*-Pi(ei)) 



(A; • p/) 

Group the terms proportional to (/* ■ Pi,/(£i)) 

= (*• • Pi(eo) fi ■ (p/(ei) - ho - a(6)] + * 



(k-Pi 
} = 



(K£)-«(£i)]-p/(£i)) V 



(k-p f ) 



I, 



+ 
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+(<*•««■)) (l ■ {»«.) - Wf)-«K6)] + t ([a(t) -ff- ft(6)) }) + 

Group in an analogous manner the terms, proportional to • •••) : 

4 (*<eo ■ {„<&) - wo - teoi + + 

4 • {»(&) - WO - + _ 1+ 

The residual terms are oc (5a) 2 : 

+{[(pm -Pf(Ci)) - 1] (/• • m ■ k) - (r • m ■ Pi (Ci))(k ■ Pf ) - (i* ■ k)(i -Pf&m -Pi)} x 

y - «(eo) 2 _ 

2(k- Pi )(k-p f ) 

Use Eq. (1301) and introduce wherever reasonable the functions of £: 

= (/* • (* • p/(0) + (** • p/(6)) • Pi(0) + 

,/,* / t «/, , x [ Q (0 - Q (£i)f , n* ft\\n m t a (0 - a (£i)] 2 , 



2(fc-P/) v ' ,vsyM 7 2(fc- Pi ) 
+ {[(p*(£i) -P/(£i)) - 1] (/* • k)(l ■ k) - (I* ■ k)(l ■ ■ pf) - (I* ■ k)(l ■ p/(fc))(* • pt)} x 

K0-«(6)) 2 , 



X 2(k- Pl )(k- Pf ) ' 

The total of all terms proportional to (a(£) — a(£i)) 2 given above vanishes. This can be 
shown using the following properties: the tensor I £g> I* is Hermitian; the difference pf — p% 
is the total of two terms, proportional to k and k', the latter being orthogonal to I. The 
residual terms are: 

#^M0 + *I^ffl) + 

The last term, —1, we transform using Eq.f TTgj) 

-\ ho - «(6)] 2 - ■ ft(eo) - i(p/(o • p/(eo) = 
= (r -^(6)) (/ -p/(0) + («• -P/(ei)) (/ -Pi(O) - ([p/(0 -Pi(0] • [p/(eo -Pi(a)]) + 

l2 /// \ /r \\2 



[g(0-q(ei)r((fc-Pt)- (fc • p/))' 

4 (A; -pt)(A; -p/) 
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On substituting this into the integral expression for the matrix element, some terms give 
zero contributions to the integral. Particularly, the following integrals vanish, as long as the 
expressions in the square brackets are the perfect time derivatives: 



{k •p i)/ (0)exp(/ , d£ 2 ) 



(k'Pu 



d£ = Q. 



From here it is also easy to derive that: 



(PM) - Pf(0 - k') exp( f {k \\ Pi ^ ] d&) 

J \ K ' Pf,i) 



d£ = 0. 



Therefore, with the transformed product of the Dirac matrices, which should be multiplied 
by a factor of two (as long as we should sum over the final electron polarization instead of 
averaging), we obtain: 



dW 



fi 



ajtft T (0 T (-^ Dd ^ 
d(k ■ k')d 2 k' ± (27rX c ) 2 (k ■ k')(k ■ Pi )(k ■ p f ) 



(32) 



where 



= (**■ •»(£)) + 



- a (£i)] 2 ((k ■ Pi) — (k ■ p/)f 



S(k ■ pi){k ■ p f ) 

The matrix element may also be summed, if desired, for two possible directions of the 
polarization vector. The second term in the integrand is simply multiplied by two, while 
in the first one the negative of the metric tensor should be substituted for the product 
of the polarization vectors (see Section 8 in [13]). so that — ■ Pi(£i)) substitutes for 

J2i (I* • Pi{€i)){l " Pi(Q)- The latter may be transformed using Eq.f fT8|) . thus, giving: 



ED 



-1 + 



[«(0-«(£i)] 2 ((fc-^ 2 + (fc-p/) 2 ) 

A{k-Pi){k-pf) 



G. Calculation of the matrix element: integration 

Now perform integration over k^. On developing the dot-product, (k' - pi), in T(£) in the 
TST metric, G^, one can find: 



T(£m-£i)=exp \i{T 1 + T 2 )} 



where 



{k-pi 



f (k-k') 
2X c (k-k')(k- Pf ) \(k- Pi ) 
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To 



(A;-fcO{(e-6) + 4[a(e 2 )-(a}] 2 cZ6} 



2X c (k- Pi )(k- Pf ) 



(a) 



ij aj6 
^-6 



(P±.) = P±oi - (a) . 



Integration over d \d ± is performed using the following formula: 



exp (iTi) d 2 \t± = n exp (iTi) d 



'(Jfc- A;') 



The right hand side is proportional to exp(zT 1 )| k± |^ 00 — 1, where the rapidly oscillating at 
large |kj_| term results in a vanishing contribution to the integral over <i£<i£i and should be 
neglected. The integrated over probability is: 

d(Jfc-Jfc') 2vrAc7(A;-^) 2 
In strong fields the following estimates may be applied: 

da 



(k-k')~\ c (k- Pi y 



da 



Now the bounds for £ + — £_ can be consistently introduced: 



\dafdt\~ 1 < - < min (V 1 (da/c^l" 1 , 1 



(33) 



Under these bounds, first, the time interval (129|) is much greater than the formation time. 
Therefore, once the double integral over d^d^i is transformed into the integral over d(£ — 
£i)d(£ + £i)/2, the span to integrate over 6?$ = d(£ — £1) can be extended towards ±00. 
Hence, the emission probability becomes linear in £ + — because only the integration over 
d(£ + £i)/2 is performed from £_ till £ + : 

dw> i (e_,e+) = wde)(e+-e-)- 

On the other hand the difference, — is to be small enough, so that the probability of 
emission within the time interval of Eq.( l29p is much less (or at least less) than unity: 

dW fl 



-dk' ± d(k' ■ k) < 1. 



(34) 



dk' ± d(k' ■ k) 

Therefore, perturbation theory is applicable. In addition, the emission probability can be 
expressed in terms of the local electric field. Note, that consistency in ( 133]) is ensured in 
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relativistically strong electromagnetic fields as long as a <C 1, with no restriction on the 
magnitude of the electromagnetic field experienced by an electron. 

Under the condition (I33p in the integral over d9 the small differences in the vector po- 
tential may be linearized: 

-[«(£) - «(£i)] 2 = [a(£) - a(£i)] 2 = 6\ 



To 



2X c (k ■ Pi )(k • p f ) 
After this the integral, 



da\ g2 (fc • Pi) 2 + (fc ■ p f f 
dCJ 4(k- Pi )(k-p f ) 



dW 



a 



d^d(k-k') 2TxX c (k-pi) 2 
by means of a substitution, z = \da./d£\9/2, may be expressed in terms of the MacDonald 
functions, using the following equations: 



+00 


"3 




/cos 


-r 






2 




-00 







2 



^ = -7^^1/3(0 



(see Eq. (8.433) in (l8|) 



2 sin 



3 z A 
2 r [J + Z , 



dz = ^K 2/3 (r), 



+00 



1 

- sin 

z 



3 z 3 
2 r [j + Z J 



2 

dz = --j= J K 1/3 (r')dr' 



as well as the recurrent relationships between the MacDonald functions: 



dKJr) 



-tf„_i(r) - K u+1 (r), K v _ x {r) - K v+1 {r) = —K v {r), 



dr r 
and the identity, K_ u {r) = K u (r). In this way we arrive at Eq.Q. The advantage of the 
MacDonald functions is the simplicity and fast convergence of the integral representation 
for them: 



KJr) 



K v {r')dr' 



exp [— rcosh(z)] cosh.(vz)dz. 



cosh(z/2;) 

exp — rcosh(z) dz, 

cosh(z) 



(see Eq. (8.432) in 18| ). In numerical simulations, therefore, the MacDonald functions are 
very easy to use. 
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VII. APPENDIX B. PHOTON IN THE QED-STRONG FIELD: ABSORPTION 
PROBABILITY 



Here we consider a photon with the wave four-vector k' and find a probability of its 
absorption in the strong ID wave field with producing electron and positron, their four- 
momenta being p e and p p correspondingly. 

The matrix element for the absorption probability is now related to the element of the 
electron phase volume. The wave function for the absorbed photon is the complex conjugated 
wave function of the emitted photon. The bi-spinor amplitude for positron is the Dirac 
conjugated amplitude for electron. With these minor changes, the matrix element for the 
absorption probability is transformed as follows: 



rfW/i (£_,£+) _ « /|_ + r(QM(p e )(7-/)M( Pp )^ 
d(k ■ p e )d 2 p e± " (47rA c ) 2 (fc ■ k'){k ■ p e ){k ■ p p ) 

T(0 = exp 



where 



and 



>>c(k-Pp,, 



(k'-Pe, P (0) 



(35) 



Eq. (l35p may be obtained from Eq. fl3TJ|) as well as the new phase multiplier may be obtained 
from the earlier derived phase multiplier by virtue of the following transformation: 



k' ->• —k', 



(36) 



Pf^Pe, Pi^-Pp, 

because the photon emission and the photon absorption are two cross-invariant channels of 
the same reaction. Then, to calculate the matrix element, one can re- write it as the double 
integral over d£d£i. The integrand in this double integral includes the phase multiplier and 
the following matrix product: 



fi(Pe(0)(7 • O«fo(0) [fi(Pe(£l))(7 • 0«(P P (6))]* 



(37) 



Now we should follow the way we used to expand Eq. ([3~Tj) with the following modifications. 
The polarization matrix for the positron is l/2[(7 • p p (£)) — !]■ While applying Eq.f l25j) to 
a positron, the dimensionless vector potential, which had been defined above in terms of a 
negative charge of electron, should now be taken with the opposite sign, therefore: 

1 



-Sp{ 



1 + S (7 - [a(ei) - a{e)]) 



x 
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x[( T -Pe(£)) + l] (7- Ox 



[(7"P P (6))-1] (7 -HI- 



We see that again the transformation Eq. (|3"6T) allows us to obtain the above matrix product 
from that one we derived in developing Eq. (1311) in Appendix A. However, in addition to the 
modifications listed in Eq. (1361) we need to transform I — > I* and change the whole sign of the 
matrix product. Now we can apply this transformation procedure to the resulting expression 
for the emission probability, Eq. (l3"2l) . and obtain the result for the absorption probability: 

d(k ■ p e )d 2 p e± (27rXc) 2 (k-k')(k-p e )(k-p p ) ' 1 j 

where 

Note an interesting polarization property: in the field of linearly polarized ID wave 
the emission probability is maximal for photons with the same polarization as that for 
ID wave, but the absorption probability for such photons is minimal (the polarization- 
dependent term in D is negative). We still employ here the probabilities averaged over the 
photon polarizations, but we think this approach may be accurate only when applied to the 
processes in circularly polarized strong wave fields. 

On integrating the averaged absorption probability over <ip e j_, we obtain the following 
expression: 

dWfi _ a (- f r °° K 5/3 {y)dy + ^g^K 2/3 {r) 



d(k ■ k>)d£ V37r\ c {k ■ k') 2 

da 



(39) 



(k ■ k') _ 3 

Xe (k-p p y Xe 2^' Pe 



d£ 



X 



c, 



which is used as the kernel in the collision integral in the present paper. 
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